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Abstract 

04 , In this article, we study the mass spectrum of the scalar and axial- vector heavy di- 

^ ' quark states with the QCD sum rules in a systematic way Once the reasonable values 

are obtained, we can take them as basic parameters and study the new charmonium- 
^ | ' like states as the tetraquark states. 
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The scattering amplitude for one-gluon exchange in an SU (N c ) gauge theory is propor- 
tional to 



^ ■ „„ „„ N c + 1„ r , , , . iV c -l 



T ki T ij = — ^j-j—($jk$ii - S ik Sji) + 4jV (SjkSu + SikSji) , (1) 



-^j- . where the T a is the generator of the gauge group, and the i,j and k, I are the color 

indexes of the two quarks in the incoming and outgoing channels respectively. For N c = 3, 
the negative sign in front of the antisymmetric antitriplet indicates the interaction is 
attractive, while the positive sign in front of the symmetric sextet indicates the interaction 
is repulsive pQ . The attractive interaction favors the formation of the diquark states in the 
color antitriplet, and the most stable diquark states maybe exist in the color antitriplet 
3 C , flavor antitriplet 3/ and spin singlet l s channels due to Fermi-Dirac statistics [2]. On 
the other hand, the study based on the random instanton liquid model indicates that the 
instanton induced quark-quark interactions are weakly repulsive in the vector and axial- 
vector channels, strongly repulsive in the pseudoscalar channel, and strongly attractive in 
the scalar and tensor channels [3]. 

The conception of diquarks has many phenomenological applications [H [6] , for exam- 
ple, the quark-diquark bound states scenario of the ground state baryons [7], the AI = \ 
rule in the non-leptonic weak decays [8], the 4 : 1 ratio of the proton and neutron deep 
inelastic structure functions in the limit x — > 1 [9], the color superconductivity in the cold 
dense quark matters [TU] . 

Taking the diquarks as basic constituents, we can obtain a new spectroscopy for the 
mesons and baryons I12|. The numerous candidates with J PC = ++ below 2GeV 
cannot be accommodated in one qq nonet, some are supposed to be glueballs, molecular 
states and tetraquark states [5], [T3l [T^] - The a(980) and /(980) are good candidates for 
the KK molecular states [15], however, their cousins cr(600) and ft(800) lie considerably 
higher than the corresponding thresholds, it is difficult to identify them as the 7T7t and ttK 
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molecular states respectively. There maybe different dynamics dominating the ++ mesons 
below and above 1 GeV which result in two scalar nonets below 1.7 GeV. The strong 
attractions between the diquark states (qq)-^ and (qq)s in S*-wave may result in a nonet 
tetraquark states manifested below 1 GeV. The conventional 3 Po QQ nonet would have 
masses about (1.2 — 1.6) GeV, and the well established 3 Pi and 3 i-2 qq nonets with J PC = 
1 ++ and 2 ++ respectively lie in the same region. Furthermore, there are enough candidates 
for the 3 P qq nonet mesons, a (1450), / (1370), if* (1430), / (1500) and / (1710) [3 H3J 

E]. 

In the tetraquark models, the structures of the nonet scalar mesons in the ideal mixing 
limit can be symbolically written as 

cr (600) = udud, / (980) 

a Q (980) = dsns, ag(980) 

K + (800) = udds, k°(800) 

The four light isospin-^ Kir resonances near 800 MeV, known as the k(800) mesons, have 
not been firmly established yet, there are still controversy about their existence due to 
the large width and nearby Kir threshold [16] . The E791 collaboration observed a low- 
mass scalar Kir resonance with the Breit-Wigner mass (797 ± 19 ± 43) MeV and width 
(410 ± 43 ± 87) MeV in the decay D + -> K~ir + ir^ [17\, and the BES collaboration ob- 
served a clear low mass enhancement in the invariant Kir mass distribution in the de- 
cay J/ip -> K*(892)K + ir' with the Breit-Wigner mass (878 ± 23+^) MeV and width 
(499 ± 52ip)MeV [32]. Recently, the BES collaboration reported the charged k(800) in 
the decay J/ip -> ^(892)^^^ with the Breit-Wigner mass (826 ± 49±||) MeV and 
width (449 ± 156+gf ) MeV [19]. 

In general, we may expect constructing the tetraquark currents and studying the nonet 
scalar mesons below 1 GeV as the tetraquark states with the QCD sum rules approach 
[201 121] . which is powerful tool in studying the ground state hadrons. For the conven- 
tional mesons and baryons, the "single-pole + continuum states" model works well in 
representing the phenomenological spectral densities, the continuum states are usually 
approximated by the contributions from the asymptotic quarks and gluons, the Borel win- 
dows are rather large and reliable QCD sum rules can be obtained. However, for the light 
flavor tetraquark states (and pentaquark states, for example, the 0(1540)), we cannot 
obtain a Borel window to satisfy the two criteria (pole dominance and convergence of the 
operator product expansion) of the QCD sum rules [22] . For the heavy tetraquark states 
and molecular states, the two criteria can be satisfied, but the Borel windows are rather 
small [231 [21|. 

We can take the colored diquarks as point particles and describe them as the scalar, 
pseudoscalar, vector, axial-vector and tensor fields respectively to overcome the embar- 
rassment [25]. In Ref. [26] . we construct the color singlet tetraquark currents with the 
scalar diquark fields, take the diquark masses as the basic parameters, parameterize the 
nonperturbative effects with the new vacuum condensates besides the gluon condensate, 
and perform the standard procedure of the QCD sum rules to study the nonet scalar 
mesons below 1 GeV, the numerical results are satisfactory. 



usus + dsds 

usus — dsds , . . - 
— y= , Oq (980) = usds, 

udus, k°(800) = usud, k~(800) = dsud. 
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In recent years, the Babar, Belle, CLEO, DO, CDF and FOCUS collaborations have 
discovered (or confirmed) a large number of charmonium-like states, and revitalized the 
interest in the spectroscopy of the charmonium states [27]. For example, Maiani et 
al identify the X(3872) observed by the Belle collaboration [28] as a 1 ++ state with 
the diquark-antidiquark structure [cq]s=i[cq]s=o + [cq] s=o [cq] s=i [2H]- The Z(4430) and 
Z(4050), Z(4250) observed in the decay modes tp'ir + and Xci^ + respectively by the Belle 
collaboration are the most interesting subjects [3D1 EH E2] . They can't be pure cc states 
due to the positive charge, and may be ccud tetraquark states (irrespective of the molecule 
type and the diquark-antidiquark type) . 

In Refs.[33j 134]. the mass spectrum of the scalar light diquark states are studied using 
the QCD sum rules. So it is interesting to study the scalar and axial-vector heavy diquark 
states with the QCD sum rules. Once reasonable values of the heavy diquark masses are 
obtained, we can take them as basic parameters and study the new charmonium-like states 
as the tetraquark states. 

There have been several theoretical approaches to deal with the heavy diquark masses, 
such as the relativistic quark model based on a quasipotential approach in QCD [35], the 
Bethe-Salpeter equation [36J, the constituent diquark model [29l ETJ [38] , etc. 

The article is arranged as follows: we derive the QCD sum rules for the scalar and 
axial-vector heavy diquark states in Sect. 2; in Sect. 3, we present the numerical results and 
discussions; and Sect. 4 is reserved for our conclusions. 



2 The scalar and axial-vector heavy diquark states with 
QCD Sum Rules 

In the following, we write down the interpolating currents for the scalar and axial-vector 
heavy diquark states, 

J\x) = e^ k qj(x)C l5 Q k (x), 
rf( x ) = e^ k sj(x)C l5 Q k (x), 
Jj(x) = e^ k qj(x)C^Q k (x), 

rf^x) = e^ k sj(x)C^Q k (x), (2) 

the i, j, k are color indexes, Q = b, c, and the C is the charge conjugation matrix. 
The two-point correlation functions H(p) and U^ u {p) can be written as 

n W = ,/A^(0|TW^(0)]|0>, 

iv (?) = ij^xe^mj^juom, (3) 

where the currents J(x) and J^ix) denote the J l (x), rf{x) and J^(x), rf„(x) respectively. 

The one-gluon exchange results in strong attractions in the color antitriplet channel 3 C , 
the quark-quark system maybe form quasibound states (diquark states) which are charac- 
terized by the correlation length L. At the distance I > L, the 3 C diquark state combines 
with the one quark or one 3 C antidiquark to form a baryon state or a tetraquark state, 



3 



while at the distance I < L, the 3 C diquark states dissociate into asymptotic quarks and 
gluons gradually. We carry out the operator product expansion at not so deep Euclidean 
space where the approximation of the correlation functions by the perturbative contribu- 
tions and the vacuum condensates makes sense. If we take the diquark state as an effective 
colored hadron and the diquark mass as an effective quantity, Mj$ ~ j-, the correlation 
function can be continued to the physical region, where the quark-quark correlations exist, 
although the diquarks are not asymptotic states, there are significant differences between 
the diquark states and conventional hadrons. The correlation functions are approximated 
by a pole term plus a perturbative continuum. 

We can insert a complete set of intermediate "hadronic" states (effective hadron states) 
with the same quantum numbers as the current operators into the correlation functions to 
obtain the "hadronic" representation [20]. Isolating the ground state contributions from 
the pole terms of the scalar and axial-vector heavy diquarks, we get the results, 

ds Ps{s) 



n(p) 


poo 
J A 2 




roo 
= / * 


Ps(s) 


= A|J(. 


Pa(s) 


= AV( 
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s — p z 



ds^4(-^ + P -^) + -, (4) 
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2 Imll^(s) 



7T 



| per 6(s - s ) , (5) 



where the pole residues Xs and (to be more precise, the current-diquark coupling 
strengths) are defined as 

<0|jM0)|S»'(p)> = XsSij, 

(0|J/^(0)|^») = Xas^, (6) 

the is the polarization vector, the A 2 is the threshold parameter, the so is the continuum 
threshold parameter, and the per denotes the perturbative contributions in the operator 
product expansion. 

In the following, we briefly outline the operator product expansion for the correlation 
functions in perturbative QCD. The calculations are performed at large space-like mo- 
mentum region p 2 -C 0, which corresponds to small distance x ~ required by validity of 
operator product expansion. We write down the "full" propagators Sij(x) and Sq(x) of a 
massive quark in the presence of the vacuum condensates firstly |21| . 



+ l -^rn s {sg s oGs) £ - —^g s G%^ + + • • • 
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where {sg s aGs) = {sg s a a/ 3G a/3 s) and ( asGG ) = { aaGa ^ G — ^ then contract the quark fields 
in the correlation functions with Wick theorem, and obtain the result: 



Substitute the full s, c and b quark propagators into above correlation functions and 
complete the integral in coordinate space, then integrate over the variable k, we can 
obtain the correlation functions at the level of quark-gluon degrees of freedom. Once the 
analytical expressions are obtained, then we can take the dualities below the thresholds 
so (the continuum states above the thresholds sq are asymptotic quarks and the spectral 
densities on both sides coincide) and perform the Borel transform with respect to the 
variable P 2 = —p 2 , finally we obtain the following sum rules for the heavy diquark states 
contain one s quark, 



where the subscript (or superscript) k denotes the scalar and axial-vector channels, i.e. 



n^(p) 



n(p) 




^{ lb S^\x) lb CSl kl {x)c} , 
{l^\x) lv CSl k ,{x)c) . 



(8) 




(9) 
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Ps(s) 



dx [x(l — x)(3s — 2mg) + m s mg] — 2rriQ(ss)5(s — tuq) 




n Q \ ~2 
—2 \S(s - m Q ) 




(10) 



5 



3 f 1 

Pa{s) = ^2 J dx [x(l - x)(2s - mg) + m s mQ] - 2mQ(ss)5(s - rag) 
m s m 2 Q (ss) m 3 Q (sg s aGs) 




fill 



ttt-q = — jp-, q = — ^, A| = (mg + m s ) 2 , and the M 2 is the Borel parameter. In Eqs.(9-ll), 
we use the dispersion relation to write the spectral densities in a compact form, the inte- 

tjfc s ~ gk s „ 

grals of the types J A ° dsf(s)e m^S(s — rriq) and J A ° dsf(s)e m^S(s — ffig) should be car- 

gk _ s _ _ m Q nS^ S 

ried out formally, i.e. J A ° dsf(s)e 'm 1 5(s — mq) = /(mg)e m 2 " and J A a dsf(s)e 'm' i 5(s — 

k k 

mq) = f{mQ)e~^ 1 despite the values A 2 , > mg and A 2 , > or < m^, where the f(s) 
denotes the spectral densities concerning the vacuum condensates. With a simple replace- 
ment m s — > 0, (ss) — > (qq), (sg s aGs) — > (qg s aGq), we can obtain the corresponding sum 
rules for the heavy diquark states contain one q quark. 

Differentiate Eq.(9) with respect to then eliminate the pole residues we can 
obtain the sum rules for the diquark masses, 

M 2 k = " J ; • (12) 

fAdsp k (s)e 

k 



3 Numerical Results 

The input parameters are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2) (qq), (qg s aGq) = m 2 (qq), (sg s aGs) = m 2 (ss), m 2 = (0.8 ± 0.2) GeV 2 , 
m s = (0.14 ± 0.01) GeV, m c = (1.35 ± 0.10) GeV and m b = (4.7 ± 0.1) GeV at the energy 
scale fi = 1 GeV [201 EH [39]. 

The Q-quark masses appearing in the perturbative terms are usually taken to be 
the pole masses in the QCD sum rules, while the choice of the thq in the leading-order 
coefficients of the higher-dimensional terms is arbitrary [40 \ [41] . The MS mass m c (m 2 ) 

relates with the pole mass m c through the relation m c (m 2 ) = rh c 

In this article, we take the approximation m c (m 2 ) ~ m c without the a s corrections for 
consistency. The value listed in the Review of Particle Physics is m c (m 2 ) = 1.27^o'iI GeV 
|16j . it is reasonable to take rh c = m c (1 GeV 2 ) = (1.35 ± 0.10) GeV. For the b quark, the 
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MS mass is mb(m 2 ) = 4.20^q; 07 GeV [IS], the gap between the energy scale p = 4.2 GeV 
and 1 GeV is rather large, the approximation mj, ~ m&(m?) ~ mft(lGeV 2 ) seems rather 
crude. It would be better to understand the quark masses m c and m& we take at the 
energy scale p 2 = 1 GeV 2 as the effective quark masses (or just the mass parameters). 
Our previous works on the mass spectrum of the heavy and doubly heavy baryon states 
indicate such parameters can lead to satisfactory results [42] , 

In the conventional QCD sum rules |20[|21j . there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter M 2 and 
threshold parameter sq. In practice, we usually consult the experimental data in choosing 
those parameters. 

Here we take a short digression to illustrate the two criteria of the QCD sum rules. 
The pole contributions (or the ratios of the pole contributions) Rk are defined by 

Q 2 dse'^pki-s) 

Rk = poo , 7T ' 

J A 2 dse m 7 Pk{s) 

k 

for a definite channel k at the hadronic representation, and the pole dominance requires 
Rk > 50%, i.e. the pole contributions dominate over the continuum contributions. At the 
level of quark-gluon degrees of freedom, the convergence of the operator product expansion 
requires the operators of increasing dimension of mass (for example, 1, 3, 4, 5, • • • ) should 
have smaller contributions. 

If we multiply the Bethe-Salpeter amplitudes of the diquark states by a charge conjunc- 
tion matrix C, the scalar and axial-vector diquark states have the same Bethe-Salpeter 
equation as the pseudoscalar and vector mesons respectively, except for the interacting 
kernels have an additional factor \ [43], the scalar and axial- vector diquark states maybe 
have slightly larger (or equal) masses than (or of) that of the corresponding pseudoscalar 
and vector mesons respectively. In the QCD sum rules for the conventional mesons and 
baryons, we usually take the energy gap between the ground states and the first radial 
excited states to be 0.5 GeV. We take the approximation Ms ~ Mp and Ma ~ My, and 
determine the central values of the threshold parameters tentatively, s s = (Mp +0.5) 2 GeV 
and = {My + 0.5) 2 GeV, where the S and A denote the scalar and axial- vector diquark 
states respectively, the P and V denote the corresponding pseudoscalar and vector mesons 
respectively. 

In calculation, we take analogous pole contributions and uniform Borel windows, i.e. 
Mmax ~ ^min = 1-0 GeV 2 and 1.5 GeV 2 in the charmed and bottom channels respectively, 
where the platforms are rather flat. The revelent parameters are shown explicitly in 
Table 1, from the Table, we can see that the pole contributions are in the range (45 — 
86)%, the contributions from the different terms in the operator product expansion have 
the hierarchy: perturbative-term 3> (qq) + (qg s aGq) 3> ( ° s *f ), the two criteria (pole 
dominance and convergence of the operator product expansion) are well satisfied. Taking 
into account all uncertainties of the relevant parameters, we can obtain the values of the 
masses and pole resides of the scalar and axial- vector heavy diquark states, which are 
shown in Tables 2-3. From Table 2, we can see that the scalar and axial-vector diquark 
states have almost degenerate masses with the corresponding pseudoscalar and vector 
mesons respectively. We should bear in mind that those values are not necessarily the 
lowest masses. 
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In this article, we intend to obtain the possible lowest masses, which correspond to 
the largest correlation lengths, and impose the two criteria of the QCD sum rules on 
the scalar and axial-vector heavy diquark states to choose the Borel parameter M 2 and 
threshold parameter so, i.e. we take smaller threshold parameters and Borel parameters 
(also Borel windows) than that presented in Table 1, and adjust them to warrant the 
uniform pole contributions (about (46 — 71)%, the smallest pole contribution presented 
in Table 1). The preferred values are shown in Fig.l and Table 4, where we can see 
explicitly that the pole contributions are about (45 — 70)%, and the contributions from 
the different terms in the operator product expansion have the hierarchy: perturbative- 
term > (qq) + (qg s aGq) S> ), the two criteria of the QCD sum rules are well satisfied 

also. On the other hand, the values of the masses and pole residues are rather stable with 
variations of the Borel parameters in the Borel windows. In fact, we can take even smaller 
threshold parameters than that presented in Table 4, however, the Borel windows are too 
small to make reliable predictions. 

Taking into account all uncertainties of the relevant parameters, finally we obtain the 
(lowest) values of the masses and pole resides of the scalar and axial-vector heavy diquark 
states, which are shown in Figs. 2-3 and Tables 2-3. In Table 2, we also present the masses 
from the relativistic quark model based on a quasipotential approach in QCD [35J, the 
Bethe-Salpeter equation [36J, and the constituent diquark model [29l ETJ [38] . From the 
table, we can see that the values from different theoretical approaches differ from each 
other greatly, and one should be careful when using them. In Ref . |26j . we introduce new 
QCD sum rules to study the nonet scalar mesons and take the values of the scalar diquark 
masses from the QCD sum rules for consistency [34J. 

The SU (3) breaking effects for the masses of the scalar and axial- vector heavy diquark 
states are buried in the uncertainties. Naively, we expect the axial-vector heavy diquark 
states have larger masses than the corresponding scalar heavy diquark states. From Table 
2, we can see that it is not the case, they have degenerate masses. Lattice QCD calculations 
for the light flavors indicate that the strong attraction in the scalar diquark channels favors 
the formation of good diquarks, the weaker attraction in the axial-vector diquark channels 
maybe form bad diquarks, the energy gap between the axial- vector and scalar diquarks is 
about | of the A-nucleon mass splitting, i.e. ~ 0.2 GeV |44[ I45j . which is expected from 
the hypersplitting color-spin interaction c Ti-TjCTi-aj, where the C is a coefficient El[5]. 
The coupled rainbow Dyson- Schwinger equation and ladder Bethe-Salpeter equation also 
indicate such an energy hierarchy [46]. Comparing with the light diquark states, the 
contribution from the hypersplitting color-spin interaction — ^— T{ ■ Tjffi ■ <jj to the heavy 
diquark states is greatly suppressed due to the large constituent quark masses, and the 
scalar and axial-vector heavy diquark states have almost degenerate masses. 

4 Conclusion 

In this article, we study the mass spectrum of the scalar and axial-vector heavy diquark 
states with the QCD sum rules in a systematic way. The diquark masses are basic pa- 
rameters in studying the tetraquark states, once reasonable values are obtained, we can 
study the new charmonium-like states as the tetraquark states with the new QCD sum 
rules developed in our previous work. 
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M 2 (GeV 2 ) 


s (GeV 2 ) 


pole 


perturbative 


(qq) + (qg s crGq) 




cg(0+) 


1.2 - 2.2 


5.6 ±0.3 


(45 - 86)% 


(59 - 73)% 


(26 - 40)% 


< i% 


cg(l+) 


1.7-2.7 


6.3 ±0.3 


(45 - 78)% 


(68 - 78)% 


(23 - 33)% 


< i% 


cs(0 + ) 


1.4 - 2.4 


6.1 ±0.3 


(46 - 83)% 


(76 - 85)% 


(15 - 23)% 


< i% 


cs(l+) 


2.0 - 3.0 


6.8 ±0.3 


(46 - 74)% 


(82 - 87)% 


(13 - 19)% 


< i% 


69(0+) 


3.8 - 5.3 


33.0 ± 1.0 


(46 - 76)% 


(64 - 71)% 


(29 - 36)% 


< i% 


69(1+) 


4.5 - 6.0 


34.0 ± 1.0 


(46 - 72)% 


(68 - 74)% 


(26 - 32)% 


< i% 


bs(0+) 


4.5 - 6.0 


34.5 ± 1.0 


(46 - 72)% 


(80 - 84)% 


(16 - 20)% 


< i% 


bs(l + ) 


4.9 - 6.4 


35.0 ± 1.0 


(46 - 71)% 


(81 - 85)% 


(16 - 20)% 


< i% 



Table 1: The Borel parameters M 2 and threshold parameters sq for the heavy diquark 
states. The "pole" stands for the contribution from the pole term to the spectral density. 
The "perturbative" stands for the contribution from the perturbative term in the operator 
product expansion, etc, where the central value of the threshold parameter so is taken. 
The contributions from the (qq) and (qg s aGq) are not shown independently for simplicity, 
and (qq) > (qg s aGq). 





M 


M 


Ref.[35j 


Refs.[29, 37, 38J 


Ref.[36 


Ref.[16j 


cq(0 + ) 


1.86 ±0.10 


1.77 ±0.08 


1.793 


1.933 


2.088 


1.867 [D] 


cq(l + ) 


1.96 ±0.10 


1.76 ±0.08 


2.036 




2.067 


2.009 [£>*] 


cs(0+) 


1.98 ±0.10 


1.84 ±0.08 


2.091 


1.955 


2.192 


1.969 [D s ] 


cs(l+) 


2.08 ± 0.09 


1.84 ±0.08 


2.158 




2.168 


2.112 [D*] 


bq(0 + ) 


5.23 ± 0.09 


5.14 ±0.12 


5.359 


5.267 


5.556 


5.279 [B] 


6,(1+) 


5.28 ± 0.09 


5.13 ±0.11 


5.381 




5.539 


5.325 [B*] 


bs(0 + ) 


5.35 ± 0.09 


5.20 ±0.07 


5.462 




5.648 


5.366 [B s ] 


6s(l+) 


5.38 ± 0.09 


5.20 ±0.08 


5.482 




5.636 


5.415 [B*] 



Table 2: The masses M of the scalar and axial- vector heavy diquark states, the values are 
in unit of GeV, the wide-hat denotes the values from the parameters presented in Table 
1. 
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A 


A 


cq(0+) 


0.53 ±0.08 


0.43 ± 0.05 


cq(l+) 


0.57 ±0.08 


0.40 ± 0.04 


cs(0 + ) 


0.64 ±0.09 


0.48 ± 0.05 


ca(l+) 


0.69 ± 0.08 


0.45 ± 0.05 


6g(0+) 


1.02 ±0.14 


0.78 ±0.14 


bq(l+) 


1.08 ±0.14 


0.76 ±0.12 


6s(0+) 


1.27 ±0.17 


0.91 ±0.14 


6a(l+) 


1.28 ±0.17 


0.88 ±0.14 



Table 3: The pole residues A of the scalar and axial-vector heavy diquark states, the 
values are in unit of GeV 2 , the wide-hat denotes the values from the parameters presented 
in Table 1. 





M 2 (GeV 2 ) 


*o(GeV 2 ) 


pole 


perturbative 


(qq) + (qg s aGq) 




cq(0+) 


1.2 - 1.8 


4.8 ±0.2 


(46 - 76)% 


(54 - 62)% 


(37 - 45)% 


< 1% 


cg(l+) 


1.3-1.9 


4.8 ±0.2 


(47 - 75)% 


(53 - 60)% 


(42 - 50)% 


<3% 


cs(0+) 


1.2 - 1.9 


5.0 ±0.2 


(45 - 77)% 


(69 - 76)% 


(23 - 30)% 


< 1% 


cs(l + ) 


1.3 - 2.0 


5.0 ±0.2 


(46 - 76)% 


(68 - 74)% 


(27 - 34)% 


< 2% 


bq(0+) 


3.5 - 4.2 


30.0 ± 1.0 


(45 - 67)% 


(53 - 55)% 


(45 - 47)% 


< 1% 


bq(l+) 


3.6-4.3 


30.0 ± 1.0 


(46 - 67)% 


(52 - 54)% 


(46 - 49)% 


< 1% 


bs(0+) 


3.5-4.4 


31.0 ± 1.0 


(46 - 71)% 


(71 - 73)% 


(27 - 29)% 


< 1% 


bs(l+) 


3.6-4.6 


31.0 ± 1.0 


(45 - 71)% 


(70 - 72)% 


(28 - 31)% 


< 1% 



Table 4: The preferred Borel parameters M 2 and threshold parameters so f° r the heavy 
diquark states. The "pole" stands for the contribution from the pole term to the spectral 
density. The "perturbative" stands for the contribution from the perturbative term in the 
operator product expansion, etc, where the central value of the threshold parameter so 
is taken. The contributions from the (qq) and (qg s crGq) are not shown independently for 
simplicity, and (qq) S> (qg s aGq). 
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Figure 1: The contributions of the pole terms to the spectral densities with variation 
of the Borel parameter M 2 . The A, B, C and D denote the channels cq, cs, bq and bs 
respectively. The /3 corresponds to the central values of the threshold parameters, the 
energy gaps among a, f3 and 7 are 0.2 GeV and 1.0 GeV for the charmed and bottom 
diquark states respectively. The (I) and (II) denote the scalar and axial-vector diquark 
states respectively. ^ 
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Figure 2: The values of the diquark masses with variation of the Borel parameter M 2 . 
The A, B, C and D denote the channels cq, cs, bq and bs respectively. The a, j3 and 7 
denote the upper bound, the central value and the lower bound respectively. The (I) and 
(II) denote the scalar and axial-vector diquark states respectively. 
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Figure 3: The values of the pole residues with variation of the Borel parameter M 2 . The 
A, B, C and D denote the channels cq, cs, bq and bs respectively. The a, (3 and 7 denote 
the upper bound, the central value and the lower bound respectively. The (I) and (II) 
denote the scalar and axial-vector diquark states respectively. 
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